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Abstract 

The two neutrino double beta decay of 94 ' 96 Zr, 98 ' 100 Mo, 104 Ru and 110 Pd nuclei for 0+ -> 0+ 
transition is studied in the PHFB model in conjunction with the summation method. In the first 
step, the reliability of the intrinsic wave functions has been established by obtaining an overall 
agreement between a number of theoretically calculated spectroscopic properties and the available 
experimental data for 94,96^94,96,98,100^ 98,ioo,i04 RU; i04,no pd and no Cd isot opes. Subse- 
quently, the PHFB wave functions of the above mentioned nuclei are employed to calculate the 
nuclear transition matrix elements M% v as well as half- lives Ty 2 . Further, we have studied the 
effects of deformation on the M2 V - 



PACS numbers: 23.40.Hc, 21.60.Jz, 23.20.-g, 27.60.+j 



I. INTRODUCTION 

The implications of present studies about nuclear double beta ((3(3) decay are far 

reaching in nature. The two neutrino double beta (2v (3(3) decay is a second order process 
of weak interaction and conserves the lepton number exactly. Hence, it is allowed in the 
standard model of electroweak unification (SM). The half-life of 2v (3(3 decay, which is a 
product of accurately known phase space factor and appropriate nuclear transition matrix 
element (NTME) M 2u , has been already measured for about ten nuclei out of 35 possible 
candidates. So, the values of NTME M 2v can be extracted directly. Consequently, the 
validity of different models employed for nuclear structure calculations can be tested by 
calculating the M 2v . The neutrinoless double beta (Ou (3(3) decay is a convenient tool to test 
physics beyond the SM. The experimental as well as theoretical aspect of nuclear (3(3 decay 
have been widely reviewed over the past years [3-19]. 

Klapdor and his group have recently reported that the Ov (3(3 decay has been observed in 
76 Ge. The results are controversial but it is expected that the issue will be settled soon j^J. 
The aim of all the present experimental activities is to observe the Ou (3(3 decay. As the Ov 
(3(3 decay has not been observed so far, the nuclear models predict half-lives assuming certain 
value for the neutrino mass or conversely extract various parameters from the observed limits 
on half-lives of the Ou (3(3 decay. The reliability of predictions can be judged a priori only 
from the success of a nuclear model in explaining various observed physical properties of 
nuclei. The common practice is to calculate the M 2v to start with and compare with the 
experimentally observed value as the two decay modes involve the same set of initial and 
final nuclear wave functions. 

In 2v (3(3 decay, the total angular momentum of four S'-wave leptons can be 0, 1 or 2, and 
is equal to the total angular momentum transferred between the parent and daughter nuclei. 
The lowest 1 + state in the final nucleus of any (3(3 decay candidate lies much higher in energy 
than the first excited 2 + state. Hence, the + transition is much less probable than the 
+ — >0 + and + — >2 + transitions. Since, the + ^2 + transition has not been detected up to 
now, the present theoretical predictions can only be checked against the + — >0 + transition 
of 2v (3(3 decay. 

In all cases of the 2v (3(3 decay for + transition, it is observed that the NTMEs 

M 2v are quenched, i.e. they are smaller than those predicted for pure quaisparticle transi- 
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tions. The main objective of all nuclear structure calculations is to understand the physical 
mechanism responsible for the suppression of the Mi v - Over the past few years, the Mi v has 
been calculated mainly in three types of models, namely the shell model and its variants, 
the quasiparticle random phase approximation (QRPA) and its extensions and the alterna- 
tive models. In the recent past, the details about these models -their advantages as well as 
shortcomings- have been discussed by Suhonen et al. [3] and Faessler et al. [13I ]. 

The shell model attempts to solve the nuclear many-body problem as exactly as possible. 
Hence, it is the best choice for the calculation of the NTMEs. However, most of the (3(3 decay 
emitters are medium or heavy mass nuclei for which the number of basis states increases 
quite drastically. Few years back, it was not possible to perform a reliable shell-model 
calculation beyond the p/-shell. Hence, Haxton and Stephenson, Jr. and Vergados ?J 
have studied the (3(3 decay of 76 Ge, 82 Se and 128 > 130 Te nuclei in weak coupling limit. Recent 
large scale shell model calculations are more promising in nature 0, 13] • The calculations 
by Caurier et al. are more realistic in which the Mi v of 82 Se is calculated exactly and those 
of 76 Ge and 136 Xe are dealt in a nearly exact manner The conventional shell model and 
Monte-Carlo shell model (MCSM) [23j have been tested against each other in case of 48 Ca 
and 76 Ge and the agreement is interestingly good. Hence, it is expected that the MCSM 
could be a good alternative to conventional shell model calculations in near future. 

Vogel and Zirnbauer were the first to provide an explanation of the observed suppression 
of M.2 V in the QRPA model by a proper inclusion of ground state correlations through the 
proton-neutron p-p interaction in the S=l, T=0 channel and the calculated half-lives are 
in close agreement with all the experimental data |24j . The QRPA frequently overestimates 
the ground state correlations as a result of an increase in the strength of attractive proton- 
neutron interaction leading to the collapse of QRPA solutions. The physical value of this 
force is usually close to the point at which the QRPA solutions collapse. To cure the 
strong suppression of M2 V , several extensions of QRPA have been proposed. The most 
important proposals are inclusion of proton-neutron pairing, renormalized QRPA, higher 
QRPA, multiple commutator method (MCM) and particle number projection. However, 
none of the above methods is free from ambiguities 13J. Alternative models, as the operator 
expansion method (OEM), the broken SU(4) symmetry, two vacua RPA, the pseudo SU(3) 
and the single state dominance hypothesis (SSDH) have their own problems 12 1. 



The basic aim of nuclear many body theory is to describe as much observed properties 
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of nuclei as possible in a coherent frame. The (3(3 decay can be studied in the same frame- 
work of many other nuclear properties and decays. Over the past years, a vast amount of 
data has been collected through experimental studies involving in-beam 7-ray spectroscopy 
concerning the level energies as well as electromagnetic properties. The availability of data 
permits a rigorous and detailed critique of the ingredients of the microscopic framework 
that seeks to provide a description of nuclear (3(3 decay. However, most of the calculations 
of (3(3 decay matrix elements performed so far do not satisfy this criterion. Our aim is to 
study the 2u (3(3 decay of 94 > 96 Zr, 98 > 100 Mo, 104 Ru and 110 Pd isotopes for 0+ -> 0+ transition 
not in isolation but together with other observed nuclear phenomena. The 2v (3(3 decay 
of 100 Mo along with the spectroscopic properties has been already studied in the Projected 



Hartree-Fock-Bogoliubov (PHFB) model using the closure approximation 25, 26]. In the 
present calculation, we have avoided the closure approximation by making use of the sum- 
mation method [27[ . Further, the HFB wave function of 100 Mo are generated with improved 
accuracy. 

The structure of nuclei in the mass region ^4^100 involving Zr, Mo, Ru, Pd and Cd 
isotopes is quite complex. With the discovery of a new region of deformation around A= 100 
by Cheifetz et al. 28|, a well developed rotational spectra was observed in several neutron 
rich Mo and Ru isotopes during a study of fission fragments of 252 Cf. The B(E2:0 + — > 2 + ) 
values were observed to be as enhanced as in the rare-earth and actinide regions. This mass 
region offered a nice example of shape transition, the sudden onset of deformation at neutron 
number iV=60. The nuclei are soft vibrators for neutron number N < 60 and quasi-rotors for 
N > 60. The nuclei with neutron number iV=60 are transitional nuclei. Thus, in this mass 
region 100 Zr, 102 Mo, 104 Ru and 106 Pd are observed to be transitional cases. In case of Cd 
isotopes, similar shape transition occurs at ^4=100. Hence, it is expected that deformation 
will play a crucial role in reproducing the properties of nuclei in this mass region ^4^100. 
Moreover, it has been already conjectured that the deformation can play a crucial role in 
case of (3(3 decay of 100 Mo and 150 Nd Q, Q. Further, all the nuclei undergoing (3(3 decay 
are even-even type, in which the pairing degrees of freedom play an important role. Hence, 
it is desirable to have a model which incorporates the pairing and deformation degrees of 
freedom on equal footing in its formalism. For this purpose, the PHFB model is one of the 
most natural choices. 

Over the past twenty years, extensive studies of shape transition vis-a-vis electromagnetic 
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properties of Zr and Mo isotopes have been successfully carried out in the PHFB model 
using the pairing plus quadrupole-quadrupole (PPQQ) interaction 3^. The success of the 
PHFB model in explaining the observed experimental trends in the mass region ^4^100 
motivated us to apply the HFB wave functions to study the nuclear 2v (3(3 decay of 100 Mo^ 
100 Ru for + — >0 + transition. Further, the success of the PHFB model in conjunction 
with the PPQQ interaction in explaining the yrast spectra, reduced transition probabilities 
B(E2:0 + -> 2 + ), static quadrupole moments Q(2 + ), ^-factors g{2 + ) of 100 Mo and 100 Ru 
nuclei as well as the Ty 2 (0 + — >0 + ) of 100 Mo j3| has prompted us to apply the PHFB model 
to study the 2u (3(3 decay of some nuclei namely 94 ' 96 Zr, 98 - 100 Mo, 104 Ru and 110 Pd for + ^0+ 
transition in the mass range 94<yl<110. 

It is well known that the pairing part of the interaction (P) accounts for the sphericity 
of nucleus, whereas the quadrupole-quadrupole (QQ) interaction increases the collectivity 
in the nuclear intrinsic wave functions and makes the nucleus deformed. Hence, the PHFB 
model using the PPQQ interaction is a convenient choice to examine the explicit role of 
deformation on NTMEs M 2u . In case of 100 Mo for + transition, we have observed 

that the deformation plays an important role in reproducing a realistic M 2v Q- 

Therefore, 

we have also studied the variation of M 2v vis-a-vis the change in deformation through the 
changing strength of the QQ interaction. 

The present paper has been organized as follows. The theoretical formalism to calculate 
the half-life of 2v (3(3 decay mode has been given by Haxton and Stephenson, Jr. j^j], Doi et 
al. Q and Tomoda Q. Hence in See. II, we briefly outline steps of the above derivations 
for clarity in notations used in the present paper following Doi et al. pj. Further, we 
have presented formulae to calculate the NTME of the 2v (3(3 decay in the PHFB model 
in conjunction with the summation method. Expressions used to calculate the nuclear 
spectroscopic properties, namely yrast spectra, reduced B(E2) transition probabilities, static 
quadrupole moments and g-factors have been given by Dixit et al. |2fij |. In Sec. Ill A, 
as a test of the reliability of the wave functions, we have calculated the yrast spectra, 
reduced B(E2:0 + — > 2 + ) transition probabilities, static quadrupole moments Q{2 + ) and 
g-factors g(2 + ) of nuclei participating in the 2u (3(3 decay and compared with the available 
experimental data. Subsequently, the HFB wave functions of the above mentioned nuclei are 
employed to calculate the M 2u as well as half-lives Ty 2 in Sec. Ill B. In Sec. Ill C, the role 
of deformation on M 2v have been studied through varying strength of the QQ interaction. 
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Finally, the conclusions are given in Sec. IV. 



II. THEORETICAL FRAMEWORK 



The inverse half-life of the 2v j3j3 decay for + transition is given by 



[T 1 2 ; 2 (0 + ^0 + )]- 1 = G 2v |M 2 ,| 



(2.1) 



The integrated kinematical factor G 2 v can be calculated with good accuracy [6J and the 
NTME M 2u is given by 

(0+|| CT T+||l+)(l+||<7T+||0+> 



2v 
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E N - (Ej + E F )/2 
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Eq + E N — Ej 



where 



E = - (E! - E F ) = + m e 



(2.2) 
(2.3) 

(2.4) 



The summation over intermediate states can be completed using the summation method 
|27| and the M 2u can be written as 



M- 



2u 



En 



0] 



^ ^ ( 1) -T— mErr 



(2.5) 



where the Gamow- Teller (GT) operator T m is given by 



(2.6) 



and 



with 



A=0 



E> 



X*- m 



D x T m =[H,[H, , [H,T m ] ](^imes) 



(2.7) 



(2.8) 



Presently, we have used a Hamiltonian with PPQQ type |32| of effective two-body interac- 
tion, which is explicitly written as 



H = H sp + V(P) + Xqq V(QQ) 



(2.9) 
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where H sp denotes the single particle Hamiltonian. The pairing part of the effective two-body 
interaction V(P) is written as 

V{P) = " (l) Y.(- l Y a+30 ~ ma ~ m " a « a * a ^ ( 2 - 10 ) 

where a denotes the quantum numbers (nljm). The state a is same as a but with the sign 
of m reversed. The QQ part of the effective interaction V(QQ) is given by 

V(QQ) = " (|) E E^H*) 44 °* S (2-11) 

where 

1 /2 

92=^) ^(M) (2-12) 



The Xijg is an arbitrary parameter and the final results are obtained by setting the x qq = 1. 
The purpose of introducing \ qq is to study the role of deformation by varying the strength 
of QQ interaction. 

When the GT operator commutes with the effective two-body interaction, the Eq. (|2.8J1 
can be further simplified to 

M 2v = Y m^^Ml (2.13) 

In the case of the pseudo SU(3) model [33-35], the GT operator commutes with the two- 
body interaction and the energy denominator is a well-defined quantity without any free 
parameter. It has been evaluated exactly for 2v (5~ (3~ - Q, Q and 2u ECEC modes 
in the context of pseudo SU(3) scheme. In the present case, the model Hamiltonian is 
not isospin symmetric. Hence, the energy denominator has not the simple form shown in 
Eq. (2.13). However, the violation of isospin symmetry for the QQ part of our model 
Hamiltonian is negligible as will be evident from the parameters of the two-body interaction 
given later. Also the violation of isospin symmetry for the pairing part of the two-body 
interaction is presumably small. With these assumptions, the expression to calculate the 
NTME M 2v of 2v (3(3 decay for 0+ -> 0+ transition in the PHFB model is obtained as 
follows. 

The axially symmetric HFB intrinsic state with K=0 can be written as 

|$o) = IW(u im + v im b\ m b\ m )\0) (2.14) 
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where the creation operators b im and b ifh are given by 

C = E C ^™ai m and 4 = E^ 1 )^^^,-™ (2-15) 

o a 

Using the standard projection technique, a state with good angular momentum J is obtained 
from the HFB intrinsic state through the following relation. 



|\I> J 



MKI 



pJ 

MK 



(2J+1) 
8tt 2 



D J MK {n)R{n)\$ K )dn 



(2.16) 



where _R(f2) and D J MK (Vt) are the rotation operator and the rotation matrix respectively. 

Finally, one obtains the following expression for the NTME M.2 V of 2v (3(3 decay for 
+ — > + transition in the PHFB model using the summation method. 



where 



2v 



y (<f ||<T.<xr+T+||<r°) 

^ E + e{n n , l n , j w ) - e(n v , l v , j u ) 



\ n {Z,N) n (Z+2,N-2) 



J 1=0 



-1/2 
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x E[( 1 +^)/S^ 2 )]^ 1 (/a 
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AT-2 



Z,N 



sinOdO 



(2.17) 



n J = 



7T 

y [cfef (l + F w / Wt )] (l + F M / Ht )] V2 d^(9)sin(9)d9 
o 



and 

™(Z,Ar),(Z+2,AT-2)(#) = 



X 



1/2 



OCt II -|- -rz,AriZ+2,Ar-2 ) 

The 7r(i/) represents the proton (neutron) of nuclei involved in the 2v (3(3 decay process 
matrices Fz,n{6) and fz,N are given by 



F Z ,n(0) - E ^m^m'S^m^m'^f^rn'^^ 



(2.18) 
. The 

(2.19) 
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fz,N — / jCjja^gCijg^mp {vi rna /ui ma ) S ma - m/3 (2.20) 

i 

The required NTME M% v is calculated using the results of PHFB calculations which are 
summarized by the amplitudes (u im ,v im ) and the expansion coefficients CV/ im . In the first 
step, matrices F n,v and f w ' v are setup for the nuclei involved in the 2v j3/3 decay making 
use of 20 Gaussian quadrature points in the range (0, 7r). Finally using the Eq. (2.17), the 
required NTME can be calculated in a straightforward manner. 



III. RESULTS AND DISCUSSIONS 



The model space, single particle energies (SPE's) and two-body interactions are same as 
in our earlier calculation on 2v (3(3 decay of 100 Mo for + — >0 + transition j^]. However, we 
have included a brief discussion of them in the following for convenience. We have treated 
the doubly even nucleus 76 Sr (N=Z =38) as an inert core with the valence space spanned by 
orbits lpi/2, 2sx/2, 1^3/2, 1^5/2, 0^7/2, 0^9/2 and 0/in/ 2 for protons and neutrons. The orbit 
lpi/2 has been included in the valence space to examine the role of the Z=40 proton core 
vis-a-vis the onset of deformation in the highly neutron rich isotopes. 

The set of single particle energies (SPE's) used here are (in MeV) e(lpi/ 2 )=-0.8, 
e(0g 9/2 )=0.0, e(ld 5/2 )=5.4, e(2s 1/2 )=6.4, e(ld 3 / 2 )=7.9, e(0^ 7/2 )=8.4 and e(0h u/2 )=8.6 for 
protons and neutrons. This set of SPE's but for the e(0/tn p), which is slightly lowered, 
has been em ploy ed in a number of successful shell model |36| as well as variational model 
calculations j3J| for nuclear properties in the mass region A=100. The strengths of the 
pairing interaction is fixed through the relation G p =30/^4 MeV and G n =20/A MeV, which 
are same as used by Heestand et al. [37] to explain the experimental g{2 + ) data of some 
even-even Ge, Se, Mo, Ru, Pd, Cd and Te isotopes in Greiner's collective model (3^. For 
96 Zr, we have used G n =22/A MeV. The strengths of the like particle components of the QQ 
interaction are taken as: Xpp = Xnn — 0.0105 MeV b~ 4 , where b is oscillator parameter. 

The strength of proton-neutron (pn) component of the QQ interaction \ vn is varied so as 
to obtain the spectra of considered nuclei namely 94 ' 96 Zr, 94,96,98,ioo MO; 98,ioo,io4 RU; io4,iio pd 
and 110 Cd in optimum agreement with the experimental results. To be more specific, we 
have taken the theoretical spectra to be the optimum one if the excitation energy of the 
2 + state E 2 + is reproduced as closely as possible to the experimental value. Thus for a 
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given model space, SPE's, G p , G n and Xvvi we nave fixed x vn through the experimentally 
available energy spectra. We have given the values of \ vn in Table I. These values for the 
strength of the QQ interaction are comparable to those su gges ted by Arima on the basis of 
an empirical analysis of the effective two-body interactions |39j . All the parameters are kept 
fixed throughout the calculation. 



A. The yrast spectra and electromagnetic properties 

In Table I, we have presented yrast energies for E 2 + to E 6 + levels of all nuclei of interest. 
The agreement between the theoretically reproduced E 2 + and the experimentally observed 
E 2 + [40] is quite good. However, it is observed that in comparison to the experimental 
spectra, the theoretical spectra is more expanded. This can be corrected to some extent in 



the PHFB model in conjunction with the VAP prescription |31[. However, our aim is to 
reproduce properties of the low-lying 2 + state. Hence, we have not attempted to invoke the 
VAP prescription, which will unnecessarily complicate the calculations. 

re experimentally observed val- 



411 ] . static quadrupole moments 



In Table II we have presented the calculated as well as t 
ues of the reduced transition probabilities B(E2:0 + — > 2' 

ri 

Q{2 + ) and the gyromagnetic factors g{2 + ) 42]. We have given B(E2:0 + — > 2 + ) results 
for effective charges e e ff =0.40, 0.50 and 0.60 in columns 2 to 4 respectively. The experi- 
mentally observed values are displayed in column 5. In case of B(E2:0 + — > 2 + ), only some 
experimentally observed representative values are tabulated. It is noticed that the calculated 
values are in excellent agreement with the observed B(E2:0 + — > 2 + ) in case of 94 Zr, 94,100 Mo, 
ioo,io4 Ru and io4p d nuclei for Geff =o.60. The calculated and observed B(E2:0 + 2+) 
values are again in agreement in case of 96 Zr and 96 Mo nuclei for e e ff =0.50. The calculated 
B(E2:0 + 2+) values for e eff =0.50 differ by 0.046 and 0.004 e 2 b 2 only in case of 110 Pd 
and 110 Cd nuclei respectively from the experimental limits. The agreement between the 
theoretical and experimental B(E2:0 + — > 2 + ) values is quite good in case of 98 Mo and 98 Ru 
nuclei for e e ff =0.40. 

The theoretically calculated Q{2 + ) are tabulated in columns 6 to 8 for the same effec- 
tive charges as given above. The experimental Q{2 + ) results are given in column 9. No 
experimental Q(2 + ) result is available for 94 ' 96 Zr. It can be seen that for the same effective 
charge as used in case of B(E2:0 + — > 2 + ), the agreement between the calculated and experi- 
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mental Q{2 + ) values is quite good for 104 Ru and 110 Pd nuclei. The discrepancy between the 
calculated and experimental values are off by 0.089, 0.14 and 0.023 e b in case of 98,100 Mo 
and 100 Ru nuclei respectively. The theoretical Q{2 + ) results are quite off from the observed 
values for the rest of nuclei. 

The g{2 + ) values are calculated with gf=1.0, g^=0.0, gj= g^=0.60. No experimental 
result is available for 96 Zr and 94 ' 96 Mo. The calculated and experimentally observed g{2 + ) 
are in good agreement for 98,100 Mo, 98 Ru, 104 Pd and 110 Cd nuclei. The discrepancy between 
the theoretically calculated and experimentally observed g{2 + ) values are 0.035, 0.021 and 
0.078 nm only for 100,104 Ru and 110 Pd nuclei respectively. The theoretical g{2 + ) value of 
94 Zr is a pathological case. The calculated g(2 + ) value is 0.121 nm while the most recent 
measured value is -0.329±0.015 nm j^j]. 

^From the overall agreement between the calculated and observed electromagnetic prop- 
erties, it is clear that the PHFB wave functions of 94 > 96 Zr, 94,96,98,ioo MO; 9 8,ioo,io4 RU; io4,iio pd 

and 110 Cd nuclei generated by fixing Xpn to reproduce the yrast spectra are quite reliable. 
Hence, we proceed to calculate the NTMEs M 2v as well as half-lives TfY 2 of 94 ' 96 Zr, 98,100 Mo, 
104 Ru and 110 Pd nuclei for + ^0 + transition. 



B. Results of 2u (3(3 decay 

The phase space factors G 2v for + ^0 + transition have been given by Boehm et al. 
for g A = 1.25 [lfj. These G 2u are 2.304xl0~ 21 , 1.927xl0~ 17 , 9.709X10" 29 , 9.434xl0~ 18 , 
9.174xl0~ 21 and 3.984xl0" 19 yr" 1 for 94 ' 96 Zr, 98 < 100 Mo, 104 Ru and 110 Pd nuclei respectively. 
However, in heavy nuclei it is more justified to use the nuclear matter value of g A around 
1.0. Hence, the experimental M 2u as well as the theoretical T^ 2 are calculated for g A =1.0 
and 1.25. 

In Table III, we have compiled all the available experimental and the theoretical results 
along with our calculated M 2u and corresponding half-lives Ty 2 of 94,96 Zr, 98,100 Mo, 104 Ru 
and 110 Pd isotopes for + — >0 + transition. We have also presented the M 2u extracted 
from the experimentally observed T^ 2 in Column 5 of Table III using the given phase space 
factors. We have presented only the theoretical T^ 2 for those models for which no direct or 
indirect information about M 2v is available to us. 

The 2v (3(3 decay of 94 Zr— > 94 Mo for + ^0 + transition has been investigated exper- 
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imentally only by Arnold |44j . who reported the limit Ty 2 > 1.1x10 yr. Theoretical 
calculations have been done by employing QRPA 

Q, OEM 0, and SRQRPA [47]. The 
presently calculated half-life in PHFB model for ^=1.25 is 7.51 xlO 22 yr, which is closer 
to the value obtained in QRPA model of Staudt et al. j3| and approximately twice of the 
lower limit given by B oby k et al. [47]. On the other hand, the calculated half-life Ty 2 in 
OEM by Hirsch et al. |46| is larger than our PHFB model value for ^^=1-25 by a factor of 
22 approximately. The predicted Ty 2 in PHFB model for g A —l.Q is 1.834xl0 23 yr. 

In case of 96 Zr, all the available experimental [18,44,48-51] and theoretical results [45- 



47,52-56] along with our calculated M 2v and corresponding Ty 2 are compiled in Table III. 
In comparison to the experimental M 2u , the theoretically calculated value given by Stoica 
using SRPA(WS) 3 is to ° smalL 0n the other hand ' the calculated half-life Ty 2 in OEM 



is quite off from the observed experimental value. The M 2v calculated by Engel et al. 
using QRPA 5^ and Barabash et al. Q using QRPA (AWS) for g A =l-0 is close to the 
experimentally observed lower limit of Wieser et al. j^J. The Ty 2 calculated by Toivanen et 
al. in RQRPA (WS) and RQRPA (AWS) are 4.2xl0 19 yr and 4.4xl0 19 yr ^ respectively 
and they are quite close to the experimental value of Kawashima et al. [48]. The predicted 
half-life Ty 2 of Bobyk et al. (47 1 has a wide range and favor al 
results. On the other hand, the T ?Y 2 p redicted by Staudt et al. 



1 the available experimental 

n 

45| is in agreement with the 
experimental result of Barabash |49J] and Wieser et al. |50l| . However, the Ty 2 calculated 
in PHFB model and in SU^)^ by Rumyantsev et al. |56J] favor the experimental values of 



NEMO |44|,|49j and Wieser et al. H for g A =1.25. 

In case of 98 Mo— > 98 Ru, no experimental result for T?Y is available so far. The theoretical 

n n 2 

calculations have been carried out in QRPA |45j, OEM |46j and SRQRPA |47j. The calcu- 
lated T?Y 2 for (7^=1-25 in PHFB model is in the range given by Bobyk et al. in SRQRPA 
model |47J. In the PHFB model for ^=1.0, the predicted half-life of 2u (3/3 decay T?Y 2 is 



n 



1/2 



1.49xl0 30 yr. The predicted Tfj 2 in QRPA by Staudt et al. [45] and in OEM by Hirsch et 
al. are larger than our predicted value for ^=1.0 by approximately a factor of 2 and 4 
respectively. 

The 2v j3(3 decay of 100 Mo for + transition have been investigated by many ex- 

perimental groups [18,51,57-65] as well as theoreticians by employing different theoretical 
frameworks [29,30,45-47,52,53,56,66-68]. In comparison to the experimental M 2u , the the- 
oretically calculated value given by Stoica using SRPA(WS) j3] is too small. The 2v (3(3 
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decay rate of 100 Mo calculated by Staudt et al. |45( and Hirsch et al. using OEM 46 1 
are off from the experimental Ty 2 . For ^=1.0, the 14^ calculated by Griffith et al. 29 1 
using QRPA model favors the results of INS Baksan 57[ and LBL |6fJ due to large error 
bar in the experimental T^%. On the other hand, the M 2v predicted by Engel et al. 5^ 
and Civitarese et al. (gtJ for ou=1.0 are in agreement with the results of LBL ((3(3], LBL 



collaboration |62], UC Irvine 63 1 and ITEP+INFN 65J due to experimental error bars. The 

n n 

values of M 2v predicted in SU(4) CTT |5(J and SU3(SPH) [66j are nearly identical and close to 



the experimental result given by Vasilev et al. |57( and ITEP+INFN |65( for ^=1.25. The 
same two M 2v for ^=1.0 are in agreement with the results of UC-Irvine ELEGANTS 
V, LBL and NEMO. Further, the value of M 2u given by PHFB model, Suhonen et al. using 
QRPA(EMP) 3(} and Hirsch et al. using SU3(DEF) 



results of Elliott et al.) 
and ITEP+INFN 



58|, ELEGANTS V [59], LBL 



66 | favor the results of UC-Irvine 
I, NEMO fell, LBL collaboration 



_65| for (^=1.25. The results of SSDH [68] are in agreement with 
the experimental half-lives of LBL 0, NEMO 0, LBL collaboration Q, UC-Irvine 



and ITEP+INFN 



The Ty 2 calculated by Bobyk et al. |4j| is in agreement with all the 



experimental results due to a large range of values from (5.04-16800) x 10 yr. 

The 2v (5(5 decay of 104 Ru— ^ 104 Pd for + ^0 + transition has not been experimentally 
investigated so far. The theoretical calculations have been carried out in QRPA 0] and 



OEM 46j. The predicted Ty 2 in QRPA by Staudt et al. |45[ is approximately one-fourth 
of our PHFB model prediction for ^=1.25 while the half life predicted by Hirsch et al. in 
OEM Q is approximately 1.31 times larger. We predict a Ty 2 for 104 Ru to be 5.73x10 
yr for ^=1.0. 

The 2v (5(5 decay of 110 Pd— s> 110 Cd for + ^0 + transition has been investigated experi 
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mentally by Winter only 
Q, SRPA(WSJ_3 and SSDH 



69| long back and theoretically by employing QRPA 



ied by Winter 



67 



45|, OEM 



The (5(5 decay of 110 Pd^ 110 Cd transition was stud- 
deducing a half-life T^ 2 > 6.0 xlO 16 yr for 2v (5(5 decay mode and a total 
half-life > 6.0xl0 17 yr for all modes. The calculated Ty 2 for =1.25 in the present PHFB 
model is 1.41xl0 20 yr, which is close to those of Semenov et al. [71] 1.6xl0 20 yr and twice 
of Civitarese et al. [67] 0.7xl0 20 yr in SSDH. On the other hand, the calculated half-life by 



Stoica |3| m SRPA(WS) is 1.186xl0 21 yr for the same g^- The calculated average half-life 



by Staudt et al. 



in QRPA is 1.16xl0 19 yr and by Hirsch et al. [46] is 1.24xl0 21 yr. For 



g A =l-0, we predict a T(J 2 for 110 Pd to be 3.44xl0 20 yr. 
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It is clear from the above discussions that the validity of nuclear models presently em- 
ployed to calculate the NTMEs M 2u as well as half-lives Ty 2 cannot be uniquely established 
due to large error bars in experimental results as well as uncertainty in Further work 
is necessary both in the experimental as well as theoretical front to judge the relative ap- 
plicability success and failure of various nuclear models used so far for the study of 2u (5(5 
decay processes. 

C. Deformation effect 

To understand the role of deformation on the NTME M 2v) we have investigated the 
variation of the latter with respect to the change in strength of the QQ interaction x m . In 
Fig. 1, we have displayed the dependence of M 2v on the Xqq f° r the 2is (5(5 decay of 94 ' 96 Zr. 
In case of 94 Zr, the M 2v remains almost constant as the strength of Xqq is changed from 0.00 
to 0.80. As the strength of Xqq is increased further up to 1.5, the M 2v decreases except at 
1.1, 1.3 and 1.5, where there is an increase in the value of M 2v . In case of 96 Zr, the M 2v 
remains almost constant as the Xqq is changed from 0.00 to 0.60. The M 2v decreases as the 
Xqq is changed to 1.20. As the Xqq is further varied to 1.5, the M 2v increases initially and 
remains almost constant. In case of 96 Zr^ 96 Mo, the experimental M 2v is available. It is 
interesting to observe that the M 2v gets tuned towards the realistic value as the Xqq acquires 
a physical value around 1.0. 

The dependence of M 2v on the Xqq has been displayed for the 2z/ (5(5 decay of 98,100 Mo 
in Fig. 2. In case of 98 Mo, the M 2v remains almost constant as the Xqq is varied from 0.00 
to 0.60 and then decreases, while the Xqq is changed to 1.2 except at 0.95. With further 
increase in Xqq, the M 2v increases at Xgg=l-3 and 1.4 and then decreases at x<?<?=1-5. In case 
of 100 Mo, the M 2v increases as the Xqq is varied from 0.00 to 0.80 and then decreases, while 
the Xqq is changed to 1.10 except at 0.95. There is a further increase in M 2v as the Xqq is 
changed from 1.10 to 1.30 and then decreases up to 1.5. It is interesting to observe that in 
case of 100 Mo^ 100 Ru, the M 2u also gets tuned towards the realistic value as the x qq acquires 
a physical value around 1.0. 

In Fig. 3, we have displayed the dependence of M 2v on the Xqq f° r the 2z/ (5(5 decay 
of 104 Ru and 110 Pd. The M 2v remains almost constant as the Xqq is varied from 0.00 to 
0.60 and then decreases as the Xqq is changed from 0.6 to 1.5 in case of 104 Ru and 110 Pd. 
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To summarize, we have shown that the deformations of the HFB intrinsic states play an 
important role in reproducing a realistic M 2u . 

To quantify the effect of deformation on M 2u , we define a quantity D 2u as the ratio of 
M 2u at zero deformation (xqq — 0) an d full deformation (xqq — !)■ The D 2v is given by 

n M 2v [Xqq = 0) m . 

The values of D 2u are 2.29, 3.70, 1.86, 2.33, 5.47 and 3.14 for 94 < 96 Zr, 98 < 100 Mo, 104 Ru and 
110 Pd nuclei respectively. These values of D 2v suggest that the M 2u is quenched by a factor 
of approximately 2 to 5.5 in the mass region 94< A <110 due to deformation effects. 



IV. CONCLUSIONS 



As a first step, we have tested the quality of HFB wave functions by comparing the theo- 
retically calculated results for a number of spectroscopic properties of 94 ' 96 Zr, 94 ' 96 ' 98 : 100 Mo, 
98 ' 100 ' 104 Ru, 104 ' 110 Pcl and 110 Cd nuclei with the available experimental data. To be more spe- 
cific, we have computed the yrast spectra, reduced B(E2:0 + — > 2 + ) transition probabilities, 
quadrupole moments Q(2 + ) and ^-factors g{2 + ). Subsequently, the reliability of the intrinsic 
wave functions has been tested by calculating the M 2v of 96 Zr and 100 Mo, for which the 2v 
(3(3 decay has already been measured. In case of 96 Zr and 100 Mo, the agreement between the 
theoretically calculated and experimentally observed M 2v as well as Ty 2 makes us confident 
to predict the half-lives Ty 2 for other nuclei undergoing 2v (3(3 decay in the mass region 
94< A <110. For 94 Zr, 98 Mo, 104 Ru and 110 Pd isotopes, the values of T*f 2 for g A =l. 25-1.00 
are (7.51-18.34) xlO 22 yr, (6.09-14.87) x 10 29 yr, (2.35-5.73) x 10 22 yr and (1.41-3.44) xlO 20 yr 
respectively. 

Further, we have shown that the deformations of the intrinsic ground states of 96 Zr, 
96 > 100 Mo and 100 Ru play a crucial role in reproducing a realistic NTME in case of 96 Zr and 
100 Mo. The NTMEs M 2u are quenched by a factor of approximately 2 to 5.5 in the mass 
region 94< A <110 due to the deformation. A reasonable agreement between the calcu- 
lated and observed spectroscopic properties of 94 ' 96 Zr, 94 > 96 > 98 > 100 ]\/[o, 98 ' 100 ' 104 Ru, 104 > 110 Pd 
and 110 Cd as well as the 2v {3f3 decay rate of 94 > 96 Zr, 98 > 100 Mo, 104 Ru and 110 Pd makes us 
confident to employ the same PHFB wave functions to study the Ou [3(3 decay, which will 
be communicated in the future. 
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Table I. Excitation energies (in MeV) of J n = 2 + , 4 + , 6 + yrast states of some nuclei in the 
mass range 94< A < 110 with fixed Gp= 30/ A, Gn= 20/ A (22/ A for 96 Zr) and e(h u/2 )=8.Q 
MeV, 



Nucleus 


X P n Theory Experiment [40] Nucleus \ vn Theory Experiment [40] 


94 Zr 


0.02519 E 2+ 


0.9182 0.9183 


94 Mo 


0.02670 E 2+ 0.8715 0.871099 




E 4+ 1.9732 1.4688 




E 4+ 1.9685 1.573726 




E 6 + 


2.7993 




E 6+ 3.3136 2.42337 


96 Zr 


0.01717 E 2 + 


1.7570 1.7507 


96 Mo 


0.02557 E 2+ 0.7779 0.778213 




E 4 + 


3.5269 3.1202 




E A+ 2.0373 1.62815 




E 6+ 


9.7261 




E 6+ 3.5776 2.44064 


98 Mo 


0.01955 E 2+ 


0.7892 0.78742 


98 Ru 


0.02763 E 2+ 0.6513 0.65241 




E 4 + 


1.9522 1.51013 




E 4+ 1.9430 1.3978 




E 6+ 


3.3098 2.3438 




E 6+ 3.6548 2.2227 


100 Mo 


0.01906 E 2+ 


0.5356 0.53555 


100 Ru 


0.01838 E 2+ 0.5395 0.53959 




E 4+ 1.4719 1.13594 




E 4+ 1.5591 1.2265 




E 6+ 2.6738 




E 6+ 2.8940 2.0777 


104 Ru 


0.02110 E 2+ 0.3580 0.35799 


104p d 


0.01486 E 2+ 0.5552 0.55579 




£4+ 1-1339 0.8885 




E 4+ 1.5729 1.32359 






2.2280 1.5563 




E 6+ 2.8790 2.2498 


no pd 


0.01417 E 2+ 


0.3737 0.3738 


110 Cd 


0.01412 E 2 + 0.6576 0.657751 




E 4+ 1.1563 0.9208 




E 4+ 1.8709 1.542412 




E e + 


2.2254 1.5739 




E 6+ 3.3865 2.479893 
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Table II. Comparison of calculated and experimentally observed reduced transition proba- 
bilities B(E2:0 + — > 2+) in e 2 b 2 , static quadrupole moments Q(2 + ) in e b and g-factors g{2 + ) 
in nuclear magneton. Here B(E2) and Q{2 + ) are calculated for effective charge e p =l+e e // 
and e n = e e ff. The g{2 + ) has been calculated for gf =1.0, g\ =0.0 and gl = g u s =0.60. 



Nucleus B{E2:$+ -> 2+) 

Theory Experiment 

e e// 

0.40 0.50 0.60 



0.40 



Q(2 + ) 
Theory 

0.50 



^(2 + ) 

Experiment^ Theory Experiment 6 



0.60 



94 Zr 


0.046 0.062 0.081 0.081±0.017 


-0.168 


-0.195 


-0.222 




0.121 


-0.329±0.01£ 






0.066±0.014 












-0.26±0.06 






0.056±0.014 












-0.05±0.05 


94 Mo 


0.148 0.188 0.232 0.230±0.040 


-0.347 


-0.391 


-0.435 


-0.13±0.08 


0.343 








0.270±0.035 


















0.290±0.044 














96 Zr 


0.044 0.060 0.078 


0.055±0.022 


-0.012 


-0.015 


-0.018 




0.254 




96 Mo 


0.265 0.335 0.413 


0.310±0.047 


-0.466 


-0.524 


-0.582 


-0.20±0.08 


0.563 








0.302±0.039 


















0.288±0.016 














98 Mo 


0.234 0.302 0.378 


0.260±0.040 


-0.439 


-0.498 


-0.557 


-0.26±0.09 


0.376 


0.34±0.18 






0.270±0.040 


















0.267±0.005 














98 Ru 


0.433 0.543 0.665 


0.411±0.035 


-0.596 


-0.667 


-0.739 


-0.20±0.09 


0.528 


0.39±0.30 






0.475±0.038 








-0.03±0.14 










0.392±0.012 















continued. 
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Table II continued. 



100 Mo 0.320 0.412 


0.515 0.511±0.009 


-0.512 


-0.581 


-0.650 


-0.42±0.09 0.477 0.34±0.18 




0.516±0.010 








-0.39±0.08 






0.470±0.024 












100 Ru 0.308 0.393 


0.488 0.494±0.006 


-0.503 


-0.568 


-0.633 


-0.54±0.07 0.355 0.42±0.03 




0.493±0.003 








-0.40±0.12 


0.47±0.06 




0.501±0.010 








-0.43±0.07 




104 Ru 0.572 0.732 


0.912 0.93±0.06 


-0.684 


-0.774 


-0.864 


-0.76±0.19 0.339 0.41±0.05 




1.04±0.16 








-0.70±0.08 






0.841±0.016 








-0.66±0.05 




104 Pd 0.361 0.460 


0.571 0.547±0.038 


-0.543 


-0.613 


-0.682 


-0.47±0.10 0.439 0.46±0.04 




0.61±0.09 










0.40±0.05 




0.535±0.035 










0.38±0.04 


110 Pd 0.479 0.614 0.766 0.780±0.120 


-0.626 


-0.708 


-0.791 


-0.72±0.14 0.47S 


0.37±0.03 




0.820±0.080 








-0.55±0.08 


0.35±0.03 




0.860±0.060 








-0.47±0.03 


0.31±0.03 


110 Cd 0.427 0.548 0.685 0.504±0.040 


-0.590 


-0.668 


-0.746 


-0.40±0.04 0.35S 


0.31±0.07 



0.467±0.019 -0.39±0.06 0.28±0.05 










45 


0±0.020 




-0.36±0.08 0.285±0.055 


"Reference 


4! 


]; b Reference 


42 


; c Reference 


43 


] 
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Table III. Experimentally observed and theoretically calculated M 2u and half-lives T^Y 2 for 
+ transition of 94,96 Zr, 98,100 Mo, 104 Ru and 110 Pd nuclei in different nuclear models. 

The numbers corresponding to (a) and (b) are calculated for =1.25 and 1.0 respectively. 



Nuclei 




Experiment 










Theory 










Ref. Projects 


^(yr) 




\M 2u \ 


Ref. Models 


\M 2u \ 




^(yr) 


94 Zr 




[44] NEMO 


> l.lxl0~ 5 


(a) 


<62.815 


* 


PHFB 


0.076 


(a) 


7.51 


(10 22 


y r ) 






(b) < 98.148 








(b) 18.34 














[47] SRQRPA 






3.08-659 














[46] OEM 






168 














[45J 


QRPA 






6.93 


96 Zr 




[50] tgch. 


0.94±0.32 


(a) 


0.074 + °-°l 7 


* 


PHFB 


0.058 


(a) 


1.56 


(10 19 


y r ) 






(b) 


llfi +0 - 027 








(b) 


3.80 






[44] NEMO 


2.1+^±0.2 


(a) 


o.o5oj£°8 


M 


SRQRPA 






0.452-61 










(b) 


0.0781°;™ 


[56] 


SU(4) CTT 


0.0678 


(a) 


1.13 






[49] NEMO 


2.0±g;§±0.5 


(a) 
(b) 


0.051+0-^ 


[55] 


RQRPA(WS) 




(b) 


2.76 
4.2 






[42] tg C h. 


3.9±0.9 


(a) 


o n3fi +0 - 005 




RQRPA(AWS) 






4.4 










(b) 


u - uo ' -0.006 


M 


QRPA(AWS) 


0.12-0.31 


(a) 


0.054-0.36 






[18] Average 


1 4+3-5 


(a) 


0.0611°;^ 








(b) 


0.13-0.88 






Value 




(b) 


0.095l° : °£ 


[53] 


SRPA(WS) 


0.022 


(a) 


10.72 






[511 Recommended 2.1+8 4 


(a) 


o.o5oj2:g8? 








(b) 


26.18 






Value 




(b) 


o.o78i8;8?^ 


[46] 
[45J 
[52] 


OEM 

QRPA 

QRPA 


0.124 


(a) 
(b) 


20.2 
1.08 
0.34 
0.82 
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table HI continued. 



98 Mo 
(10 29 yr) 



* PHFB 



[47] 
[4i 

14; 



SRQRPA 

OEM 

QRPA 



0.130 (a) 6.09 
(b) 14.87 
4.06-15.2 
61.6 
29.6 



100 Mo [65] 


ITEP+INFN 


7.2±0.9±1.8 


(a) 





i 91 +0.032 
±zi -0.018 


* 


PHFB 


0.104 


(a) 


9.79 


(10 18 yr) 






(b) 





1 qn+ - 050 

iyu -0.028 








(b) 


23.90 


M 


ITEP 


8.5 


(a) 





112 


w 


SSDH 




(a) 


7.15-8.97 








(b) 





174 


[47] 


SRQRPA 






5.04-16800 


[63] 


UC Irvine 


6.82l°;|±0.68 


(a) 
(b) 






1 9P-+0.013 
lz,J -0.009 

iOK+0.020 
ly,J -0.014 


[56] 


SU(4) (7T 


0.1606 


(a) 
(b) 


4.11 
10.03 


w 


LBL+MHC+ 


7.61 2 ; 2 


(a) 





1 1 O+0.013 
iio -0.014 


m 


SSDH 


0.18 


(a) 


3.27 




UNM+INEL 




(b) 





1 oc;+0.020 
io,J -0.022 








(b) 


7.99 


[61] 


NEMO 


9.5±0.4±0.9 


(a) 
(b) 






1 nfi +0.008 
±uu -0.007 

1fi i-+0.013 
iud -0.010 


[53] 


SRPA(WS) 


0.059 


(a) 
(b) 


30.45 
74.34 


[60] 


LBL 


9.7±4.9 


(a) 
(b) 






1 nc; +0.044 

iud -0.019 

1 £o+0.069 
±uo -0.030 


[66] 


SU(3)(SPH) 


0.152 


(a) 
(b) 


4.59 
11.2 


M 


ELEGANTS V 11.5±|:g 


(a) 





nofi+0.010 

uyo -0.011 




SU(3)(DEF) 


0.108 


(a) 


9.09 








(b) 





1I5n +o.oi5 

i,JU -0.016 








(b) 


22.19 


[58] 


UC Irvine 


1 1 6 +3 - 4 

il -°-0.8 


(a) 





nOR+0.004 
uyu -0.012 


[46] 


OEM 






35.8 








(b) 





1 40+0.005 

± ^ y -0.018 


[30] 


QRPA(EMP) 0.101 


(a) 


10.39 


[57] 


INS Baksan 


q q+2.0 
°-°-1.0 


(a) 





1 7q +0.036 
1 ' J -0.038 








(b) 


25.37 








(b) 





98n+ - 055 

^ OU -0.059 


[29J 


QRPA(EMP) 0.256 


(a) 


1.62 


[IS] 


Average 


8.0±0.6 


(a) 





1 1 P-+0.005 
iid -0.004 








(b) 


3.95 




Value 




(b) 





1 on+0-007 
iou -0.006 


[45J 


QRPA 






1.13 


[51] 


Average 


8.0±0.7 


(a) 





1 1 C-+0.005 
iid -0.005 


[52] 


QRPA 


0.211 


(a) 


2.38 




Value 




(b) 





1 cn+°- 008 
iou -0.007 








(b) 


5.81 
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Table III continued. 



104 Ru * PHFB 0.068 (a) 2.35 

(10 22 yr) (b) 5.73 

[46j OEM 3.09 

[45j QRPA 0.629 



110 Pd [69] > 6.0xl0~ 4 (a) < 6.468 * PHFB 0.133 (a) 1.41 
(10 20 yr) (b) < 10.106 (b) 3.44 

[71] SSDH 1.6 
[67] SSDH 0.19 (a) 0.7 

(b) 1.70 

[20] SRPA(WS) 0.046 (a) 11.86 

(b) 28.96 
[46] OEM 12.4 
[45] QRPA 0.116 



gch. denotes geochemical experiment and *Present work. 
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